NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 





TECHNICAL NOTE 2226 


THEORETICAL ANALYSIS OF OSCILLATIONS IN HOVERING 
OF HELICOPTER BLADES WITH INC LINED AND OFFSET 
FLAPPING AND LAGGING HINGE AXES 
By M. Morduchow and F. G. Hinchey 
Polytechnic Institute of Brooklyn 



Washington 
December 1950 


i ta i L l i 


r - ■ i 
r La 



LIBRARY KAFB, Nllfl 



TECH UBRARY KAFB, NM 

11 


0D£>5 


□1 


NATIONAL ADVISORY COMMUTER FOR AERONAUTICS 


TECHNICAL NOTE 2226 


THEORETICAL ANALYSIS OF OSCILLATIONS IN HOVERING 
OF HELICOPTER BLADES WITH INCLINED AND OFFSET 
FLAPPING AND LAGGING HINGE AXES 
By M. Morduchow and F. G. Hinchey 


SUMMARY 


The frequency and damping characteristics of the coupled flapping 
and lagging oscillations of helicopter blades in hovering are derived 
for the general case in which the lagging (vertical) hinge axis is off- 
set from the flapping hinge axis, while both hinge axes are inclined. 

An offset of the flapping hinge axis from the rotor axis of rotation is 
also considered. An inclination of the hinge axes causes pitch change 
of a blade section due to angular displacements about these axes. It 
is found from this analysis that an inclination of the flapping hinge 
axis chiefly affects the values of the natural flapping frequencies and 
the damping factors of the lagging motion of a blade. An inclination 
of the lagging hinge axis, which causes considerable flapping motion 
during the lagging oscillations, mainly affects the natural frequencies 
and damping factors in lagging. The analysis and the numerical examples 
worked out here indicate that significant increases in the damping of 
the lagging motions, which ordinarily border on instability, can be 
obtained by suitable inclinations of the hinge axes, especially of the 
lagging axis. Offsetting the flapping and lagging hinge axes tends 
especially to increase the natural lagging frequency. 


INTRODUCTION 


The aim of this investigation is to determine theoretically the 
characteristics of the rigid-body oscillations of helicopter rotor 
blades with inclined and offset flapping "(horizontal) and lagging 
(vertical) hinge axes. For this purpose a simple method which requires 
the solution of a fourth-degree algebraic equation for obtaining the 
natural frequencies and damping decrements of the small coupled flapping 
and lagging oscillations of a blade in hovering about a state of static 
equilibrium is developed here . __ Approximate formulas for these charac- 
teristics, based on the quartic equation, are, moreover, given. 
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The equations of motion of a "blade are derived hy means of the 
Lagrange equations for a system having two degrees of freedom, -namely, 
rotation about an inclined lagging hinge axis and rotation about an 
inclined flapping hinge axis. The flapping axis is assumed as offset 
from the rotor axis of rotation, while the lagging axis -is further 
ass ume d as offset from the flapping axis. The aerodynamic loads in 
hovering are derived here on the basis of quasi-stationary flow, with 
the assumption of constant downwash throughout the rotor disk. 

Basically, the effect of inclining either hinge axis is to create 
a kinema tic relation between flapping or lagging and change of pitch of 
a blade section. This kinematic relation is of the form 0 = f(P,£), 
where 0, p, and £ are, respectively, the pitch, flapping, and lagging 
angles. It may be noted that, 'although the effects of several kinematic 
constraints have already been derived in reference 1, the condition 
mentioned above has not been explicitly investigated there. Since the 
aerodynamic loads in hovering are largely functions of the blade, or 
pitch, angle, it follows that in the present analysis the effect of 
inclination of the hinge axes appears essentially in the resulting 
effect on the aerodynamic loads caused by the changes in blade angle. 

The analysis is divided into three sections. In the first section 
the explicit equations of motion of a blade are derived. The equations, 
pnd their solutions, determining the steady-state (static in hovering) 
values of flapping, lagging, and blade . angles are included here . 

In .the second section the dynamic equations are solved in general 
terms for the frequency and damping characteristics of the oscillations 
of a blade. Coupling between flapping and lagging, which becomes espe- 
cially large when the lagging hinge axis is "inclined, is taken into 
account. In addition to the development of a straightforward exact 
procedure, and explicit approximate solutions, for the natural frequen- 
cies and damping decrements of the oscillations, stability criterions 
are given. For any given data, these stability criterions, as well as 
the explicit equations for the details of the oscillation characteristics, 
can be readily applied. General conclusions regarding the effect of incli- 
nation and offset of the hinge axes are derived from the general solutions 
obtained. 

Tn the final section numerical examples are worked out to illustrate 
the general conclusions reached as well as the method of application of 
the general solutions developed here. 

This investigation has been part of a project carried out at the 
Polytechnic Institute of Brooklyn under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
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SYMBOLS 


A,D 


amplitude of flapping and lagging oscillations, 
respectively 


ao,a^,a 2 ,a 2 ,a] 1- coefficients of quartic equation ( 26 a), defined by 

equations ( 26 b) 


c do 


profile-drag coefficient of blade section 


c 

c o 

e l 


e 2 

do 

^tDX* d-QyJ d-Q2 


F 1' F 2' F 3 ,F 4'T. 

F 5 ,F 6 , F ? , p 8 , J 

g 


chord length of blade section 

value of c of innermost blade section 

distance of flapping hinge axis from axis of rotation 
flapping hinge eccentricity (fig. l) 

distance of lagging hinge axis from flapping hinge 
axis; lagging hinge eccentricity (fig. l) 

section .profile drag per unit length of blade 

drag load components in x-, y-, and z-directions, 
respectively 

dimensionless, definite integrals defined by equa- 
tions ( 21 b) 

acceleration due to gravity 


H = mo/pirc 0 Z 

i - 

L* 

l 1' l 2 ,l 3 ,l 4 ,l 5 

V>V'V 

z 

M = g/ft 2 Z 


lift loads based on quasi -stationary flow 

dimensionless definite integrals defined by equa- 
tions ( 21 b) 

lift load components per unit blade length in x-, y-, 
and z-directions, respectively 

length of blade measured from lagging hinge axis 



k 

m 

mo 

n 

P 

Ewt 

P x >Py»Pz 

Qi 

Qp>Q£ 

q 

q± , 
qf^qz 

^fo^^Zo 

^^zi 

\ 

R 

r 

s 

S 1 

T 

t 

V 
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mass per unit length of blade 

value of m at innermost blade section 

number of rotor blades 

complex frequency; if p = ft(-R ± Icd) (R and o> real), 
then 0Jfl/2jt is the natural frequency in cycles 
per second, while 2itR/oi is the logarithmic decrement 

uni t weight loads of blade in z-direction 

aerodynamic plus weight loads pqr unit blade length 
acting in x-, y-, and z -directions, respectively 

generalized force 

generalized forces corresponding to generalized 
coordinates p and £ 

dimensionless complex frequency .(p/ft) 

generalized coordinate 

dimensionless flapping and lagging complex frequencies, 
respectively 

dimensionless uncoupled flapping and lagging complex 
frequencies, respectively 

first approximation to dimensionless coupled flapping 
and lagging complex frequencies, respectively 
(equations (28) and (30)) • . 

tip radius of blade 

distance of any point on blade from axis z of 
rotation 

coordinate along blade axis measured from lagging 
hinge * 

value of s for innermost blade section 
kinetic energy of system 
'time 

vector of resultant inflow velocity 
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'ns 


resultant inflow velocity in plane of "blade section 


V x >Vy,V z 


W 

w i 

X,Y,Z 


x,y,z 

x = x + (flr) x 
y = y + (flr) y 

h 

0 

0,1 


components of resultant inflow velocity in hovering 
at any point of "blade 


gross weight of helicopter 

i 

induced downwash 

forces along axes of rotor system- 

2 n 2 


P 2 


X = 


P x dsi 


S 1 


Y = 


p y dsj Z = 


P z ds 


ls l 


rotating right-handed Cartesian coordinate system in 
which z is axis of rotation and x coincides with 
undeflected longitudinal "blade axis (figs. 1 and 2) 


magnitude of angular displacement about inclined 
flapping hinge axis 

component of 0g about a horizontal axis, (AE in fig. l) 
normal to deflected blade axis (OC) 

deviations of 0 and _ £ from steady-state position 


0o>£o>®o 


steady- state values of 0, . and 0, respectively 
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circulation vector . • 

inclination of lagging (vertical) hinge axis measured 
from a line parallel to z-axis- (see fig. 2), positive 
when positive lagging (i.e., retreating) produces a 
decrease in blade angle ' 

inclination of flapping (horizontal) hinge axis measured 
from a line parallel to y-axis (see fig. l), positive 
when positive flapping produces a decrease in blade 
angle 


e i - e i I 1 
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G g = ®2/^ 

5 

e 

e o 

(A0)i 

(A0) 2 

5 = s/Z 

*1 = B l/ Z 
P 
0 

fl 

(ar) x ,(flr) y 


Oi 


magnitude of angular displacement about inclined 
lagging hinge sods 

component of £g about an axis (OD in fig. 2) parallel 
to z-axis (see fig. .1) 

total pitch, or blade, angle of blade section after 
flapping and lagging deflections 

design, or initial, blade angle 

change of pitch due to rotation £g about lagging 
hinge axis 

change of pitch due to rotation 0g about flapping 
hinge axis 


density of air 

angle between x-axis and radius of a point on blade to 
z-axis 

angular velocity of rotor system 

components of linear rotational velocity of blade in 
x- and y-directions, respectively 

natural frequency of oscillation of rotating blade 


BASIC EQUATIONS 


The basic equations of motion of a rotating helicopter blade in 
flapping and lagging are derived here. 

Coordinate system . - A right-h a nd e d Cartesian coordinate system 
(x,y,z) in which the z-axis coincides with the axis of rotation of 
the rotor system and the x-axis is in the direction of the undeflected 
longitudinal blade axis is used. The x- and y-axes rotate with a blade 
about the z-axis. The flapping hinge axis of a blade, which is assumed 
to be in the plane of rotation (x,y), is considered as offset by a 
distance e from the axis of rotation and inclined at an angle 5g 
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to a line parallel to the y-axis (see fig. l) . The lagging hinge axis, 
which is assumed to he in the vertical (x,z) plane, is considered to 
he offset from the flapping hinge axis hy a distance eg and inclined 
at an angle 8]_ to a line parallel to the z-axis (see fig. 2). 

Lagrang e equations . - The Lagrange equations of motion can he 
expressed in general as follows: 


dt^i) c% ^ i 


( 1 ) 


where t is time, T is the kinetic energy of the system, q^ is a 
generalized coordinate, and Q-j_ is the corresponding generalized force. 
The dot indicates differentiation with respect to time. 


The kinetic energy of a moving hlade is given hy 
# 

iZ 

I 

T = - I mV 


(■^ 2 . -i2 .2^ , 

+ y + z J ds 


(2a) 


where 


x » x + (ftr) x 1 
y = y + (flr)y 


(2b) 


2 is the length of the hlade measured from the lagging hinge, s-j^ is 

the distance of the innermost hlade section from the lagging hinge, m 
is the mass per unit length of a hlade, ds is a blade-length element, 
ft is the an g ular velocity of the rotor system, r is the distance of 
any point of the hlade from the axis of rotation, and (ftr) x and (ftr) y 

denote components of the linear rotational velocity of a hlade element 
in the x- and y-directions, respectively. 

The position of a hlade after lagging and flapping can he determ i ned 
as follows: 

% 

First consider a rotation of magnitude about the lagging hinge 

axis. This produces a change of pitch (A0)j_ of magnitude 


(A0)! = -Cs sin 5], 


(3) 
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Moreover, the component £ of rotation, about an axis parallel to the 
z-axis (OD in fig. 2) has the magnitude 

£ = £5 cos 81 (^) 


Let the blade now rotate with an angular displacement of 3g about 
the inclined flapping axis. Then the additional pitch change (&9)q 

due to this rotation is the component of rotation about the deflected 
blade axis (OC in fig. l) . Thus 


(A 0 ) 2 = 35 sin(£ - 83) 


(5) 


Moreover, the campohent 3 of rotation of the blade about a horizontal 
axis (AE in fig. l) normal to the deflected blade axis (OC) after lagging 
is given by 


3 = 35 cos (83 - £) 


(6) 


The coordinates of a point- P on the longitudinal axis of a blade at 
a distance s (fig. l) from the lagging hinge axis can be expressed in 
terms of the "flapping" and "lagging" displacement-components 3 and £ 
by the following expressions: 

3 COS 8 ? 

x = e-, + s cos 3 cos £ + e 0 cos 

1 2 cos(b 3 - 0 

y = -s cob p sill £ f 


3 cos 8-3 

z = s sin 3 + e 2 sin — 

003(83 - £) 1 

The velocity components at any point of a blade which are due to 
the rotation of the rotor system can be expressed as follows (see fig. 3 ): 


(fir) x = -ftr sin 0 = -Sly 
(nr) y = Hr cos 0 = ftx 


( 8 ) 


where x and y are given by equations ( 7 ). 
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By considering 3 and £ as the two generalized coordinates, 
equations (2a), (2fb), (7)> and (8) maybe substituted into equation (l) 
to give the left side in terms of 0, £, and time derivatives thereof. 

Generalized forces .- The generalized forces Qi here are due to 
the aerodynamic and weight loads. The force (here Q*' and Qr) 

is defined by 


where 


Q£ 8£ + Q p 83 = X 8x + Y 8y + Z 8z 


X = 


Y = 


Z = 


p x ds 


Py ds 


P z 48 


(9a) 


(9b) 


and p x , Py, and p z are the aerodynamic plus weight loads per unit 
blade length in the x-, y-, and z -directions, respectively. 

From equation (9a) it follows that 


8x 8v 8z 

= X ~ + Y ~ + Z 

8p 83 83 


a „ _ y dx Y 8y 8z 

Q ^- X 5? + X sr + Z 5r 


( 10 ) 


3x 

The derivatives and so forth can he determined in terms of 0 

and £ by means of equations (7). 

Lift loads .- The lift loads,. based on quasi- stationary flow, can 
be found from the three-dimensional Kutta- J oukowski theorem: 


l 1 = pv x r 


(11) 
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where p is the density of the fluid medium, V is the vector of 
the resultant inflow velocity, and r is the circulation vector* The 
velocity components in hovering at any point of a "blade can he expressed 
as: 


V x = fly - x 

Vy = -Ax - y ■ 
V z = -wi - z 


( 12 ) 


where w i is the induced downwash at the rotor, while x, y, and z 
can he put in terms of p and £ hy means of equations (7)» The 
induced downwash is here assumed as constant throughout the rotor disk 
in hovering. According to the simple momentum theory, it can he calcu- 
lated from the expression: 


Wi _ / W 

m V 2xR 2 pfl^ 2 


(13) 


where W is the gross weight "of the helicopter and R is the tip radius 
of a hlade. 

By considering the circulation vector . F to coincide with the 
deflected longitudinal hlade axis and to have a magnitude determined hy 
the Kutta condition of finite trail ing- edge velocity in two-dimensional 
flow, the lift loads can he derived in a manner similar to that given 
in detail in reference 1. By assuming p, £, 6 (the hlade angle), 

and wi/flZ to he first order small, the following expressions for the 
lift load components per unit length of a hlade can thus he derived to 
second-order quantities : 


L x ’ = prtcft(e^ + e 2 + s)p[^P(e 2 + s) + Wi - fl(ei+ e 2 + s)pj 

Ly’ = pirc ^(e 2 + s) + wj [^P(e 2 + s) + wi - fl(ei + e 2 + s)0^ 

L z ' = pitc ||^(e 2 + s) +w i ][ts-fl(e 1 + e 2 + s)] - 2^fls(e]_+ e 2 + s)0 

fl 2 ( e i +e 2 + s) [P£(e! + e 2 ) - ( ei + e 2 + s)©]J 
where c is the chord length of a hlade section. 
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Drag loads .- The expressions for the lift load components include 

•tiie effect of induced drag, since they include the effect of downwash. 
The profile-drag loads can he determined from the following expression 
for the profile drag per unit length of a "blade: 

*o = c do 2 ^ns^ c ( 15 ) 


where c do is the profile-drag coefficient of. a blade section and V ns 

is the resultant Inflow velocity in the plane of the blade section. By 
considering c^ 0 as constant and the resultant drag force as acting 

parallel to the resultant inflow velocity, the following expressions, 
with c do considered as first order small, are obtained for the drag 
load components to second-order quantities: 


d ox = -§ CPCCd<^ 2 s(ei + e 2 + s) 

' = 2 pcc dc^( e l + e 2 + s ) [5^® ■ fl ( e l + e 2 + s )] ► 

d oz =-| P cc do*K e l + eg + s)^ + P(e 2 + s)] 

w 

Weight loads .- The unit weight loads of a blade act in the 
z-dlrection and may be expressed in the foim: 


(16) 


= ““8 dT) 

where g is the acceleration due to gravity. 

Equations of motion .- The total pitch angle of a blade section 
after deflecting is, according to relations (3) to (6), given by: 


0 = So + P tan(£ - 63) - £ tan 81 (l8) 

where 0 Q is the initial, or design, blade angLe before flapping or ' 
lagging. This expression for 0 must be substituted in equations (l4) 
for the lift loads. 

The equations of motion of a blade in flapping and lagging can now 
be written explicitly in terms of 3 and £ by substituting into 
equation (l). By putting ' 


P = 3o + P(t) 

£ = So +T(t) 
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where j3 0 and t >Q denote the values of p and £ in the steady 

state, while p(t) and £(t) denote small deviations from the steady 
state, and observing that the steady-state values must alone satisfy 
the equations of motion, the following nondimensional equations are 
obtained for the values of p Q and to and for the vibrational 
deviations f3(t) and £(t): 


Po 1 ( 5^ e 2 + £)( 6 1 + e i 



^( € 1 + c 2 + £)( e 2 + 5 ) |( 6 i + 6 2 ) Mo + 

+ 2 r) " ^ €l + e 2 + £) e oJ d & = 0 


( 19 a) 


Cn - 


J 5i + + !) 2 + m( e i + e 2 + E) 0 o - (mf 


i d! 


n i 

H ( 6 l + 6 2) . jjM d| 
Jh “° 


(19b) 


■|(Fi - HF 2 ) + ■§ F 3 + ~2 + ^( F 5 + F 8 taQ 5l ) + + Fj) =0 (20a) 


— - y 

|(HF2 + Li) + I L 2 + — HL 3 + p(tan 83 - £ 0 sec 2 & 3 )L 4 + 

£ jk.5 + ^tan 81 - P 0 sec 2 83^L4j = 0 


where 


S 1 -T 


H = 


hiq 


P«C 0 Z 


► _ 8 1 

T 

e _ e 2 

2 = T 
M = -S— 


(20b) 


(21a) 
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Fl = 


■ Jt x ^ (62 + l){ [ 2(e i + s 2 + s > e ° - 5l] 

f 

Sx ”° ! 


dl 


p s - 2 Po J ( gjS («2 + I) « 




■^( € 2 + t) 2 ( e l + e 2 + I) d| 


F 4 “ 5^ e 2 + £) 2 d l 

F 5 = Po ~( 6 1 + e 2 + 6) ( + 6) [ 6 1 + 6 2 - ( e l + «2 + I) sec^J dl 

P 6 = J gi ^(«2 + + 6 2 + 6 “ M$o) d * 

F 7 = + 6 2 + l)( 6 2 + I) [( 6 1 + e 2) So + 

(6^ + 63 + £)(tan 83 - £ 0 sec%3)J d| ^(21b) 

F 8 = -§^ «i + e 2 + |) 2 (€ 2 + I) d| 

11 = /ei ^ + E,5 [ 2 S - ( e i + 5 a + s ) 8 0 ] « 
l2 ' f ^/l 1 ^^ 6 l + e 2 + 5 ) i! 

^ / s > 2 « 

14 = S / ix^ E ( ei + £ 2 + 8 ) a 8 
L 5 - ( e l + %> « 
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dIq and c 0 refer to the mass and chord/ respectively, of the innermost 
hlade section, and 0 O is the steady-state value- of 0 . The value 
of 0 O is given hy: 

Bo = 6 0 + 3 o ta11 ( So ~ 83) - £ 0 tan &i (22) 

In deriving equations (20a) and (20h) the vibrational deviations p 
and £ (not P Q and £ 0 ) have_been assumed sufficiently small so that 
only first powers of {T and £ may he retained. : 

The condition that the lift load support the gross weight of the 
helicopter in hovering leads to the following expression fpr 0 O : 




e n = 


nl^prtc 


fe f-( e l+ e 2 +|) 
1=1 c ° « 


l 


(23) 


. %( *1 + % * if 


where n is the number of blades in the rotor system. 


Steady-state values of p, g, and 0 .-, The steady-state equa- 
tions (19a), (19b), and (23) can be solved for the values of P 6 , £ 0 , 

and 0 O as follows. Equation ( 23 ) may be solved first for 0 O with 
the higher-order term (in . Po^o) there neglected. - The value of £ 0 can 

then be obtained from equation (l9b) and the value of P D can subse- 
quently be determined from equation ( 19 &) • A second approximation for 
the value of 0 O can now be calculated by inserting the valiies qf p Q 
and £ 0 , thus obtained, into equation ( 23 ) « Second approximations 
for P 0 and £ 0 , by using the new value of 0 O , can then be determined 
from equations (19a) and (19b). It will usually be found that the 
second approximations in this iterative procedure are sufficiently 
accurate for practical purposes. The value qf the design angle 0 q can 
be found from the relation 


00 = e o “ Po tan (£0 “ s 3) + £0 ta A & 1 
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It may "be noted here that the values of 0 O , £ Q , and 0 O are all 

independent of the inclinations 8]_ and 83 of the hinge axes. However, 
the appropriate value of the design blade angle . 6q does depend on the 

values of 6^ and 83. 


Dynamic equations - physical significance .- Equations (20a) and 
(20b) are the equations determining the characteristics of the small 
vibrations of a helicopter blade about the state of static equilibrium 
given, as explained above, by equations (19&) and (19b). The physical 
significance of each of the terms in equations (20a) and (20b) is as 
follows. The terms proportional to 0 and £ represent the inertia , 
loads per unit of length acting on a blade. The term proportional to 0 
in equation (20a) represents primarily the aerodynamic damping load 
caused by a change of angle of attack induced by a vertical velocity 0 
of a blade element. The term 0HFg in equation (20a) is due to the 
centrifugal loads, while the term 0Fy represents the lift load. The 

term |jFi in equation (20a) represents a lift load component (including 
the effect of downwash) caused by a lagging velocity £, while the 
term represents a Coriolis load. The term proportional to £ 

in equation (20a) is a lift load. 


In equation (20b), the term is proportional to the profile 

drag of a blade section. The term in this equation represents 

the centrifugal restoring force in lagging, while the term proportional 
to £Lij. represents an induced drag, due to the downwash w^. The 

A A 

term represents a Coriolis load while the term jjLi is an aero- 

dynamic load. Finally, the term 0Llj. is an aerodynamic load due to 
the downwash w^. 


Equation (20a) corresponds primarily to the flapping oscillations, 
while equation (20b) corresponds primarily to the lagging oscillations. 

It will be observed, however, that these equations are coupled with each 
other. It may be further noted that all of the coupling terms in lagging 

(proportional to 0 and 0 in equation (20b)) have magnitudes which 
are at most first order small but that in flapping there, is a finite- 
order coupling term proportional to tan 81. Consequently, inclination 
of the vertical (lagging) hinge tends to produce considerable coupling 
between flapping and lagging. 
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GENEBAL SOLUTIONS FOE STABILITY AND VIBBATION CEAEACTEEISTICS 


The general solution of the flapping and lagging equations and the 
general implications of this solution are given here. In addition, 
stability conditions are shown and discussed. 


Basic solution .- Equations (20a) and (20b) are a system of linear, 
homogeneous, differential equations with constant coefficients. The 
solution of this system can be obtained by setting 


3 = AcP* 
I = DeP t 


( 2k ) 


where A and D are real or complex constants (amplitudes) and p is 
a real or complex constant ("frequency"). 

The following linear, simultaneous, algebraic equations in A 
and D are obtained by substituting equations (24) into equations (20a) 
and (20b): 

A^q%F^ + qF 3 + HFg + F^ + olq^ - HFg) + F^ + Fg tan = 0 (25a) 


A 




+ 


^tan 63 


£0 sec% 3 )Li f ] 


+ 


D Iq^HLg + qLg + HL5 + ^tan 5^ - P Q sec 2 &q^Li|j = 0 


(25b) 


where q = p/fi. 

For the existence of a nontrivial solution to the above equations 
the determinant of, the coefficients of this system must vanish. 
Expanding this dete rminant and setting it equal to zero leads to the 
following quartic equation (to second powers of small-order quantities) 
for the complex dimensionless frequency q: 


4 

&4q + 


agq-3 + agq^ + a^q 


+ a„ = 0 


(26a) 
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where 

Ej,. = 

e^ = H^FjjLg + ^ 


&2 


&1 

(HF 2 - Fj^Lij. tan 83 - (hf 2 + L l)( F 5 + Fq tan 83 
a,0 = (BF 6 + F y)jHL 5 + (tan ^ - 0 Q sec^^L^J + 
Li,.p8 tan 8 l(^ 0 s ec 2 8 3 - tan 83) - F5 tan 6 


Frequency and damping characteristics .- Equation (26a) determines 
four values of ql If these values constitute two pairs of complex 
conjugate roots, then they can be interpreted physically as follows. 

If q '= -R ± las, where R and co are real positive quantities, 
then a ifi is the natural frequency of oscillation, while 2n=§ is the 
logarithmic decrement, that is, the natural logarithm of the ratio of 
amplitude of oscillation at any time to the , amplitude one period or 
cycle later. If R is negative, then the oscillations will increase 
in amplitude with time and the motion will be unstable. A negative 
real root for q indicates a nonoscillating subsident motion, which is 
stable, while a positive real root, indicates a nonoscillating motion of 
steadily increasing amplitude, which is unstable. 

Equation (26a) can be solved exactly for the values of q by - 
means, for example, of Ferrari’ s method (reference 2, pp. 51 and 52). 
For purposes of a quicker, but approximate, calculation, the following 
method may be used. 



= HF 4 |hL 5 + (tan - 0 Q sec^L^ + HL^HFg + F^) + FgLg 

* 

(HF 2 - F !)(hF 2 + L x ) 


= f 3 [hl 5 + ( 


tan 8 -l - 0 Q sec‘ 


-S3K] + *2^*6 + F ?) 


> (26b) 
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Let a first approximation to a pair of roots of equation ( 26 a) "be 
that which would he- obtained in flapping if there is no coupling with 
lagging. If one of these roots is denoted by q-f 0 , then it follows from 
equation (25a) with D = 0 that: 


«fb = 


£i_ 

2HFlj. 



HFg + Ff 
HF4 


(27) 


A second approximation q^ to equation ( 26 a), by using Newton* s method 
(reference 2, pp. 96-97) "With q fQ as first approximation whether it is 
real or complex, will then be: 




qfo(*i - HFg) + Fg q £o (HFg + Lj) + Lfr tan BgJ +V3 tan 8i[qf 0 (HFg +I»i) + (tan 63 - £ 0 sec^Bg)^ 


(28) 


Similarly, a first approximation to the second pair of roots q^ Q 

of equation ( 26 a), obtained by considering lagging without coupling to 
flapping (equation (25b) with A = 0), maybe assumed as 

a = j£_±l HI 5 + (tan - p Q sec g 5 3 )L 4 

q 2o 2 HL 3 M\2hl 3> / HL 3 v 

A second approximation, by Newton's method, will be 

q = l [ t Uo( F l - ^2) + *5] [llofes +L l) +L k tan +r 8 tan OiCqicfHFg + L].) +(tan 83- £ p sec 2 !^!^] ^ 

^ tC ”™~~ 


Equations (27) and (29) may be interpreted physically as giving 
the oscillation characteristics in flapping and lagging, respectively, 
if there is no coupling between these motions. Equations (28) and (30) 
give a first approximation to the effect of the coupling on the damping 
and frequency characteristics of the corresponding motions. If the 
corrections indicated by equations (28) and (30) are found to be rela- 
tively large, then it is suggested that the roots of equation (26a) be 
calculated by Ferrari ' s method. 


Corresponding to any root of equation ( 26 a) , equation ( 25a) or 
equation (25b) deter min es the value of the ratio A/d, whose absolute 
value gives the ratio of the flapping amplitude to the lagging ampli- 
tude. If there is no coupling, then 1^1 = 00 for flapping, while 

I A I ^ * 

U| = 0 .for lagging. A complex value for A/D indicates a phase 
.difference between the component flapping and lagging motions. 
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From equations (27) and ( 28 ), it will be found that in ordinary 
cases the angle of inclination Si of the lagging hinge axis has rela- 
tively little effect on the damping and frequency characteristics of the 
motion of a hlade corresponding primarily to flapping. The angle of 
inclination 83 of the flapping hinge axis has, however, an important 
effect on the frequency characteristics of the primarily flapping motion 
of a blade. Increasing 83 negatively diminishes the natural frequency 
in flapping, since the total effective restoring force (proportional 
to p in equation (20a)) is diminished. It maybe noted that a nega- 
tive 83 indicates an inclination of the flapping hinge axis such that 

a positive (upward) flapping motion produces an increase in the blade 
angle 6 (cf. equation (5))* For sufficiently large negative values 
of 83, in fact, the restoring force becomes negative and an unstable 
motion results, evidenced by at least one positive real root for q. 
Physically, the reason for the decrease of restoring force in flapping 
with increasing negative values of 83 is that a negative 83 tends 

to produce an aerodynamic (lift) load component opposing the centrifugal 
restoring force. 

From equations ( 29 ) and ( 30 ) it will be found that increasing the 
angle 8j positively 1 - increases the natural frequency of the motion 
corresponding primarily to lagging but decreases the damping factor 
(i.e., tends to make the real part of q positive). The natural fre- 
quency of the primarily lagging motion is only slightly affected by the 
angle 83- However, the damping factor in lagging is considerably 

influenced by 83 (cf. numerical examples discussed in a following 
section) . 

From equations (27) and (29), as well as equations ( 28 ) and (30), 
it will be found that, in general, oscillating motions corresponding 
to flapping will be much more heavily damped than those corresponding 
to lagging, since the negative real parts of q^ will be of higher order 

of magnitude than those of qj, 

The effect of offset, expressed quantitatively by' the parameters Gj 
and is included in the general solutions given here. It will be 

found from equation (29), in fact, that the chief effect of small off- 
sets is an increase in the uncoupled lagging frequency with increase in 
the value of G-^ + 62- 

The coupling effects, as given approximately by equations ( 28 ) 
and (30) , are much more important in the lagging motions (equation (30)) 

•*-A positive value of 81 means that the lagging hinge axis is 
inclined such that a lagging (retreating) motion of the blade produces 
a decrease in the blade angle 6 (cf. equation (3)). 
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than in the flapping motions (equation (28)). This is due to the fact 
that, although the corrections due to coupling in hoth of these equa- 
tions are at most first order small, the magnitude of the uncoupled 
lagging complex frequency q^ (especially its real part) is also small; 

consequently, the correction in lagging is relatively important. The 
relatively important effect of coupling in lagging is also evidenced hy 

the relatively large value of j^j (order of magnitude of unity) corre- 
sponding to q^ when 5^ is not zero. 

General stability criterions .- If one is interested only in the 
stability of the motion of a helicopter blade, then it is not necessary 
to determine the actual values of the four roots of equation (26a). 

The general condition that the motion be stable is that all of the 
coefficients of equation (26a) have the same sign (viz, positive), while 
also (reference 3) 

a l( a 2 a 3 “ a i a h) ~ a o a 3 2 >0 (31a) 


Since a^ and a^ are necessarily positive, it is necessary for sta- 
bility that 



a x >0 


(31b) 



Conditions (31a) and (31b) together are necessary and sufficient for 
stability. 


For purposes of a quick, but rough, calculation, the stability 
criterions for the uncoupled flapping motions may be considered. From 
expression (27) for q^ 0 , this condition is found to be approximately: 
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In actual cases, this condition will usually he satisfied, for 83 > - 45°, 
approximately. Condition (32) may he considered, as an approximation even 
with coupling taken into account, since, as noted, above, the effect of 
coupl ing is relatively small in flapping. It must he observed., however, 
that inequality ( 32 ) is only a necessary hut not a sufficient condition 
for stability of the entire motion of a blade, since it does not include 
the motions corresponding primarily to lagging. 


NUMERICAL EXAMPLE 


A numerical example is given here to indicate the type of results 
which may he obtained in practice and to illustrate the method outlined 
in the previous section for finding the oscillation characteristics of 
a rotor-blade system. The inclination of the lagging hinge axis 8 ^ 
will first be permitted to assume a range of values, while the flapping 
hinge axis will not be inclined; that is, 83 = 0°. Then the flapping 

hinge axis will be inclined at various angles 83 while the lagging 
hinge axis inclination 81 is fixed at 0 . Finally, the case of 

inclining both hinge axes simultaneously will be considered. 

The following data for' a typical helicopter are assumed; 


Gross weight of helicopter, W, pounds 3000 

Number of rotor blades, n 3 

Angular speed of rotor system, ft, radians per second ........ 25 

Tip radius of a rotor, R, feet 21.5 

Length of blade proper, Z, feet 20 

Length of eccentricity, e^, feet 1 

■Length of eccentricity, e£, feet 0.5 

Chord width at root, c Q , feet 1 

Airfoil section -.. Clark Y 

Profile-drag coefficient, c^, ...... 0.01 


The structure of the blade consists of a tubular steel spar with plywood 
ribs flanged on and plywood skin attached to them. The diameter of the 
spar at the root is 0.15 foot -and the wall thickness is 3/16 inch. The 
effective area of the wood is assumed as four-tenths that of the steel, 
and the density of steel is 15.12 slugs per cubic foot while the density 
of wood is 1.5 slugs per cubic foot. It will be further assumed that 
the blade has a constant cross section. Therefore, 

c _ m 
c D mo 


1 
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It is necessary to determine the values of H, M, and. l first 

in order to obtain the steady-state values of p, t, , and 9. By substi- 
tuting into equations (13) and (21a) one obtains 

H = 0.77^014 

M = 0.002576 

St - 0 - 0ta665 

The values of p Q , So* and 9 Q can now be found by the iterative 

procedure outlined in the section "Basic Equations." The following are 
the values obtained by approximately solving (i.e. , neglecting higher- 
order terms in B 0 and £ 0 ) equations (23), (19b), and (19a), 

respectively: 


e Q = 0.122583 
S Q = 0.051955 


p 0 = 0.070968 


A second approximation for 9 0 is now obtained by substituting the 
above values of P Q and So into equation (23) and now retaining second- 
order terms. By using this new value of 9 o, second approximations 
for So and p 0 are obtained from equations (19b) and (19a), respec- 
tively. The following are the values of the second approximations: 

t 

@ 0 = 0.122969 
S 0 = 0.052162 
p Q = 0.071369 

It can be seen that the second approximations here give sufficient 
accuracy, since they differ little from the first approximations. 
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It is now possible to calculate the parameters 


given by equations (21b). Thus, it is found that 

F-l = 0.055504 

F 5 = 0.002062 - 0.022350 sec 2 8 3 

F 2 = 0.049364 

Fg = 0.385102 

F 3 = 0.294369 

Fj = 0.001507 + O.313161 tan 83 

F4 = 0.358958 

Fg = 0.313162 


in terms of 83 


- O.OI6335 sec^ 


^ = -0.006139 
Lg = O.OOO875 


L3 = 0.333333 
L4 = 0.015451 


L5 = 0.037500 


The values of the uncoupled complex frequencies in flapping and 
lagging (qf Q and qi QJ respectively) are first obtained by solving 

equations (27) and (29), respectively. Then the first approximation to 
the coupled frequencies and q^, respectively) of the coupled 

oscillations is found from equations (28) and (30). In addition, the 
quart ic equation (26a) has been solved for q by Ferrari’ s (exact) 
method to test the accuracy of Newton’s method as it is applied here for 
finding the approximate values of the coupled complex frequencies. 
Finally, the ratio of the flapping amplitude to the lagging amplitude 
A/D can be found from either equation (25a) or (25b). It should be 
mentioned here that the ratio A/D serves as a check, since a partic- 
ular value of q should yield the same result from both equations (25a) 
and (25b). It is sometimes necessary to adjust the root slightly in 
order to satisfy both equations. 


Tables I and II contain the values of q fQ , q fl , q^ (exact), 
and A/D, and qj Q , q ?1 , qj (exact), and A/D, respectively, when 8 1 
has a range of values and 83 = 0°. Tables III and IV contain the same 
information for 5p = 0° and a range of values for 63. The results 
given in tables V and VI were obtained by inclining the lagging and 
flapping hinge axes simultaneously. 


The results of the numerical examples given in tables I to VI 
illustrate the general conclusions reached in the prece ding section. 

For example, the negative real parts of the flapping complex frequencies 
are in all cases much greater than those of the lagging complex fre- 
quencies, an indication of much greater damping of the flapping motions 
than of the lagging motions. Table II illustrates how increasing 8^ 
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. algebraically adversely affects the stability of the lagging motions , 
since the real part of q^ becomes less and less negative' and' even 
positive. Increasing 8^ negatively, on the other hand, has a bene- 
ficial effect on the damping of the lagging motions. From tables I 
and II, in fact, it may be seen that the case 83 = 0°, 81 = -k 5 ° is 

quite satisfactory from a stability viewpoint, since the damping factor 
for the lagging osc illa tions is over ten times as great as when 83^ = 0°, 

while the flapping oscillations remain stable. The increase of the 
natural lagging frequency imaginary part of. q^j with algebraically 

-tnr-rgnfi-tng 5^ is also shown here. As evidenced by the values of A/D 
in table II, considerable flapping motions accompany the lagging motions, 
xinless S x £ 0°. Table I indicates the rather small effect which an 

inclination Sj of the lagging hinge axis has on the flapping 'oscil- 
lation characteristics. 

Table IH illustrates the increase of flapping frequency with alge- 
braic increase of the an gl e 83. It also shows the danger of insta- 
bility when 83 has an algebraically low value. Here, for 83 = - 1 * 5 °, 

the positive real root for qf indicates instability. Table IV shows 
the relatively little influence of the inclination 83 of the flapping 
hinge axis on the natural lagging frequencies. It is interesting to note 
here, however, that inclining the flapping hinge axis at an angle of 
63 = 130° has a beneficial effect oh the stability characteristics of 

the lagging motion, since the real part of q^ then has a negative value 
of at least three times that of q 3 corresponding to the other values 
of 83 considered here. However, it may be observed from table III 
that, in this case, the flapping frequency (imaginary part of qf) would 
then become dangerously low. 

Of the three cases considered in tables V and VI ^83^ 0°, 83 ^ 0°^ 

the case 81 = -30°, 83 = 30° appears to be the most suitable from the 

point of view of stability characteristics. The case of 81 = 30°, 

83 = -30° evidently leads to . instability, since the real part of q^ 
is positive. Moreover, the case 81 = -30°, S3 = -30° feives a slight 
instability in flapping (q^ has a small positive real value). 

A comparison in the tables between corresponding values 
°f q^. with q f , and q^ with qj, shows that in most of the 

numerical cases treated here, Newton* s method (equations (28) and (30)) 
gives fairly good second^ approximations for the complex roots of 

^The first approximations, namely the uncoupled complex frequencies, 
are given by equations (27) and (29). 
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equation ( 26 a). 
gave unreliable 
S 3 = -30°. The 
from the tables 
or q 2 . 


The only cases here in which this second approximation 
results are: Sp = 0°, 63 = -4-5° and Sp = -30°, 

effect of coupling on the complex frequencies can be seen 
by comparing q fo with q fl or q f , and q Zo with q 21 


CONCLUSIONS 


From an investigation to determine theoretically the characteristics 
of the rigid-body oscillations in hovering of helicopter blades with 
inclined and offset flapping and lagging hinge axes, the following con- 
clusions are drawn: 

1. The frequency and damping characteristics of the (in general,' 
coupled) flapping and lagging oscillations of a helicopter blade can be 
calculated directly either by solving a quartic equation for the dimen- 
sionless complex frequency or by using a second approximation for the 
roots of this quartic. 

2. From the solutions it is- found that the main effect of inclining 
negatively the flapping (horizontal) hinge axis is to diminish the natural 
flapping frequency by di min ishing the total restoring forces in flapping. 

3. Inclination of the flapping hinge axis also has an important 
effect on the damping characteristics of the lagging motion, as 
evidenced by the numerical example treated here, which indicated an 
inclination of -30° to be a desirable value from this viewpoint. This 
value, however, gave a dangerously low natural flapping frequency. 

4. The chief effect of inclining negatively the lagging (vertical) 
hinge axis is to decrease the natural lagging frequency but to increase 
the damping factor, and therefore the stability, in lagging. Inclination ' 
of the lagging hinge axis increases the coupling between lagging and 
flapping and leads to large flapping motions accompanying the 1 a g g l n g 
motions . 

5. The main effect of a slight offset of the hinge axes is an 
increase in the natural lagging frequencies. 


of Brooklyn 
October 24, 1949 


Polytechnic Institute 
Brooklyn, N. Y., 
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TABLE I— VALDES OF q fo , q a , q f (EXACT) 
WHEN Bi HAS A RANGE OF VALDES AND 83 


8 1 

(aeg) 

Ifo 

(equation. (27)) 

Ifl 

(equation (28)) 

(exact) 

a/d 

^5 

-0.5297 ± 0.85961 

-0.6388 ± 0.93231 

-0.5858 ± 0.90381 

5.181 T 6.21*01 

30 

-.5297 ± 0.8596i 

-.5677 ± 0.88271 

-.5630 ± 0. 88161 

l*.90l* T 6.1871 

0 

-.5297 ± 0.85961 

-.5338 ± 0.84841 

-.5255 ± 0.85151 

1*.1*70 T 6.3271 

-30 

—5297 ± O.85961 

-.1*878 ± 0.821*91 

-.1*795 ± 0.821*11 

3 V 888 T 6.2171 

-^5 

-.5297 ± 0.85961 

-.1*999 ± 0.83111 

-.1*392 ± 0.80911 

3-514 T 6.11*31 


TABLE H.- VALDES QF q 2o , q u , q z (EXACT), AND A/D 
WHEN 8n HAS A RANGE OF VALDES AND 60 = 0 ° 


8 1 

(deg) 

«Zo 

(equation (29)) 

q 21 

(equation (30)) 

(exact) 

a/d 

45 

-O.OOI69 ±0.41001 

0.05671 ±0.37531 

0.05435 ± 0.38451 

-0.9307 ± 0.42281 

30 

-.00169 ± 0.37791 

.03117 ± 0.36381 

0 

.03151 ± 0.36601 

-.5346 ±0.21401 

0 

-.OOI69 ±0.32901 

-.0051*91 ±0.33181 

-.005891 ±0.33l6i 

.06312 T 0.047131 

-30 

-.00169 ±0. 27141 

-.04303 ±0.25761 

-.04303 ±0.25761 

.7447 T 0.24251 

-1*5 

-.OOI69 ±0.21991 

-.07733 ±0.20521 

-.07733 ± 0.20511 

1.302 T 0.28961 


TABLE 


IH— VALDES OF q fQ , q^, q f (EXACT), AND A/D 
WHEN 83 HAS A RANGE OF VALDES AND 82 = 0 ° 


8 3 

(dfig) 

Ifo 

(equation (27)) 

q.jfi 

(equation (28)) 

If 

(exact) 

a/d 

45 

r0.5297 ± 1.3441 

-0.5271 ± 1.3391 

-0.5271 ± 1-3391 

7.881 T 9.0421 

30 

—5297 ± 1.1701 

-5273 ± 1.1651 

-.5274 ± 1.1651 

6.444 T 8.2871 

0 

-5297 ± 0.85961 

-.5338 ± 0.84841 

-.5255 ± 0.85151 

4.470 T 6.1281 

-30 

-.5297 ± 0.26171 

-.4695 ± 0.19741 

-.5093 ± 0.18091 

3.478 T 1.0331 

-45 

f .1389 

•.1817 

.1737 

3.024 


l -1.198 

-1.222 

-1.221 

7.302 
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TABLE IV.- VALDES OF q jQ , q n , %l (EXACT), AMD A/D 
WHEN 53 HAS A RANGE OF VALDES AMD 8j_ = 0 ° 


8 3 

(cleg) 

*Io 

(equation. (29)) 

(equation (30)) 

*1 

(exact) 

a/d 

45 

-O.OOI696 ± 0.32241 

-0.004358 ± 0.32991 

-0.004360 ± 0.32981 

0.07376 T 0.023361 

30 

-.OOI696 ± 0.32681 

-.004068 ± 0.33121 

-.004165 ± 0.33111 

.05907 T 0.026801 

0 

-.OOI696 ± 0.32901 

-.005491 ± 0.33181 

. -.005891 ± 0.33161 

.06312 T 0.047131 

-30 

-.OOI696 ± 0.32681 

-.02232 ± 0.33451 

-.02202 ± 0.33701 

.08841 T 0.24081 

-*•5 

-.OOI696 ± 0.32241 

-.04357 ± 0.32241 

-.007629 ± 0.35971 

.2364 T 0.21861 


TABLE V.- VALDES OF q fo , q a , q f . (EXACT), AMD A/D 
WHEN Bx ^ 0 °- AND 83 A 0 ° 


8l 

83 

Ifo 



If 

a/d 

(cleg) 

(deg) 

(equation (27)) 

(equation (28)) 

(exact) 

30 

-30 

-0.5297 ± 0.26171 

-0.5471 ± 0.55121 

-0.6048 ± 0.40481 

4.093 T 2.4631 

-30 

30 

-.5297 ± 1-1701 

-.5073 ± 1.1441 

-.5055 ± 1-1491 

6.206 ? 8.3721 



■ 



[-8179 

1-997 

-30 

\ 

-30 

-.5297 ± 0.26171 

-.3374 ± 0.012681 

< 

1 

[ .01907 

5.260 


t 


TABLE VI.- VALDES OF q jQ , q u , q z (EXACT), AND A/D 
WHEN 8x ^ 0 ° AND 83 ± 0 ° 


8 1 

(deg) 

8 3 

(de|) 

4lo 

(equation (29)) 

*11 

(equation (30)) 

(exact) 

a/d 

30 

-30 

-0.001 696 ± 0.37601 

0.07318 ± 0.34081 

0.07341 ± 0.35891 

-0.6374 

± 0.83611 

-30 

30 

-.OOIG96 ±0.26871 

-.02806 ± 0.29891 

-.02598 ± 0.29951 

.4698 

± 0.10461 

-30 

-30 

-.OOI696 ±0.26871 

-.2173 ± 0.18201 

-.1320 ± 0.45221 

.04182 T 2.1061 














































